Upper bounds in phase synchronous weak coherent chaotic attractors 
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An approach is presented for coupled chaotic systems with weak coherent motion, from which it is 
estimated the upper bound value for the absolute phase difference in phase synchronous states. This 
approach shows that synchronicity in phase implies synchronicity in the time of events, characteristic 
explored to derive an equation to detect phase synchronization, based on the absolute difference 
between the time of these events. 



INTRODUCTION 

This work deals with the phenomenon of phase syn- 
chronization (PS) [11,3 in coupled chaotic systems, which 
describes interacting systems that have a bounded phase 
difference, despite of the fact that their amplitudes may 
be uncorrelated. PS was found in many natural and 
physical systems [HjB , being experimentally observed in 
electronic circuits 0, in electrochemical oscillators [3], 
in the Chua's circuit [6| , and in spatio-temporal systems 
0. There are also evidences of PS in communication 
processes in the Human brain p?, ^ . 

In the case of two coupled systems, PS exists [l[ if 



il < 



(1) 



where 0i.2 are the phases calculated from a projection of 
the attractor onto appropriate subspaces in which 

the trajectory has coherent properties P, [lo]. The ra- 
tional constant g [3] is the frequency ratio between the 
average phase growing and g is a finite constant to be 
determined, bounded away from zero. 

The purpose of this work is to give an upper bound 
value for the absolute phase difference in Eq. (H]) in phase 
synchronous states, in terms of a defined phase [I4I . This 
is equivalent to determine an inferior bound value for the 
constant g. We show that this minimal value, namely (r), 
can be estimated as the average growing of the phase, cal- 
culated for typical trajectories, in one of the subspaces. 
Particularly, (r)= (W) x (T), where (W) is the average 
angular frequency associated to a subspace X, and (T) 
is the average returning time of trajectories in this same 
subspace, calculated from the recurrence of events of the 
chaotic trajectory. Similarly to periodic oscillating sys- 
tems, in which it is valid to say that an angular frequency 
oj is related to the period T by w = 2tt/T, for chaotic sys- 
tems it is valid to say that (W) = {r)/{T). 

In the derivation of the constant (r) , we obtain a series 
of inequalities that can be used to check for the existence 
of PS. A particular interesting one is suitable for systems 
where the only available information is a series of time 
events. We also introduce the phase of a chaotic attractor 
to be given by the amount of rotation of the tangent 
vector of the flow. 

These results are shown to be valid to weak coherent 



attractors. By weak coherent attractors we mean follow- 
ing Ref. [2^, attractors in which it is possible to define 
a Poincare section or a threshold that defines an event, 
such that for the time between two events r, it is true 
that |t — (r)| < K, where (T) is the average returning 
time between two successive events, and k < (T) is a 
small constant. So, our results are extended to attractors 
whose trajectories might not have a clear rotation point, 
but still presenting a weak coherent property in the time 
between events, e.g. bursting/spiking dynamics. 

For illustrating our ideas, we use two coupled Rossler 
oscillators, and two coupled neuron models from the 
Rulkov map This last example was chosen because 
we want to demonstrate that PS can be detected by only 
knowing the time at which bursts occur (events). 



A MINIMAL BOUND FOR THE CONSTANT g 

We start by developing some ideas to give a minimum 
bound for c in Eq. ([1]). For simplicity, we eliminate the 
rational constant q given by q— , by a changing of 
variable, 02 (i)' = 902 (i)- With a slight abuse of notation, 
from now on, we omit the / symbol in the phase. Note 
however that such a changing of variable does not change 
the fact that PS exists or not. 

Having two oscillators 5*1 and S2 that are coupled form- 
ing the attractor E, we define the susbspaces Xj to be a 
special projection in the variables of S. This projection 
is such that the attractor in this subspace presents the 
coherent properties defined in ji^l . Subspaces Xj are the 
same ones where the phase is calculated. 

Next, we define a time series of events, where events 
here are the crossing of the trajectory to a given Poincare 
section, some local maxima/minima, or the crossing of 
one variable to some threshold. Being E the attractor of 
the coupled system, and Xi and X2 two subspaces (on 
which the phase is defined), r* is the time at which the 
i-th event happens in Xj . We consider the average return 
time, (Tj), of the subspace Xj to be the average of time 



intervals Tj=Tj 
events. So, 



between two events in Xj, for N 
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We introduce the phase as the amount of rotation of 
the unitary tangent vector, Aj{t). Being \Aj{t + St) — 
Aj{t)\ a small displacement of the phase for the time 
interval St, calculated on the subspaces Xj, and making 
St — > 0, we arrive at 



\A,\dt 



(3) 



So, (j>j{t) measures how much the tangent vector of the 
flow, projected on the subspaces Xj, rotates in time. This 

equation also suggest that \Aj\ can be seen as an angu- 
lar frequency, more precisely Wj — 
average angular frequency is simply 



dt 



Ml 



and the 



1 



W^dt. 



(4) 



We introduce the quantity r\ — j^f Widt, which is 

the evolution of the phase from the time T2 (when the i-th 
event happens in X2) until the time T2^^ (when the (i+l)- 
th event happens in X2). Thus, (ri) = (X^i^o ''1)/-^' 
in a continuous form, after the iV-th event, this average 

i + l 

is calculated as (ri) — J^i Wjdt/N which is equal 
to 



{ri) = 



lo' Widt 



N 



(5) 



Using that it is valid to say that ~ N{Ti). For r, 
00, in Eq. ^ we have that (ri) = 77^^ /g^ Wi 

("^ lo^ M^idi^ , which using Eq. (O, can be written 



dt 



iTT) 

as 



(Ti] 



(6) 



These calculations can be done for (r2), however, if 
PS exists, i.e. Eq. ([T]) is satisfied, one should have that 
{W,MW2), (ri) = (r2), and {T,) = (T^). Thus, in Eq. 
we can use the index j. 

Synchronicity of events: The number of events at 
a given time for synchronous oscillators is not always the 
same, but can differ by an unity. This occurs because 
the A^-th event in Xi and X2 may not be simultaneous, 
resulting in a difference of an unity between the number 
A^i and N2 of events, in Xi and X2, respectively. So, 
we can say that the number of events in PS are always 
related by 



\Ni{t) ~ N2{t)\ < 1. 



(7) 



The inequality in Eq. ([7|) is another variant for an 
equation already used to detect phase synchronization 
In that equation, every time an event occur, like the 



crossing of the trajectory through a treshold, the phase 
is assumed to grow 2tt. And PS is considered to happen 
if the phase difference is always smaller or equal than 27r. 
Note that Eq. ([7]) can also be used to detect synchronous 
events in maps, in the case an event can be well specified. 
As an example, one can observe the occurrence of local 
maxima in the trajectory [l6j . 

Synchronicity in the time of events: 

Using Eq. ([2]) in Eq. ([7]), we arrive at: 



N 



(8) 



This equation is related to the weak coherence in the 
dynamics. The more phase coherent the attractors are 
the more the amount | J^ii^^i ~ ^2)! approaches to zero. 
As a consequence, the value (Ti) over estimate the max- 
imum difference in the time intervals between events. 
To overcome this, we introduce a physical parameter, 
namely 7, which brings us information about the coher- 
ence of a specific system. Thus, we put Eq. (|5]) as 



iV 
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1=0 



(9) 



It is important to notice that 7 also brings some infor- 
mation about the projection and about the section in 
which the events are defined, once that the difference in 
the time intervals depends on the projection and on the 
Poincare section definition. Our calculations to the cou- 
pled Rossler-like attractors, show that 7 — 1/2. 

Multiplying both sides of Eq. ^ by (W^i), we can 
relate time of events with the averaging growing of the 
phase: 



N 



{w^)\j2{Ti-n)\<j{n). 



(10) 



i=0 



Synchronicity of the phase: Next, we represent 
Eq. (HI) at the time the A^-th event happens in Xi, by 



N-l 

i=0 



<{r), 



where 



Widt 



W2dt 



(11) 



(12) 



The term ^fjQ^{r\ — r\) represents the phase in Xi at 
the moment the (A^— l)-th event happens in X2 minus the 
phase in X2 at the moment the {N — l)-th event happens 
in Xi. The term S,{N) represents the difference between 
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the evolution of the phase from the event A'' — 1 till the 
time at which the 7V-th event happens at the subspace 
Xi, minus the evolution of the phase at the subspace 
X2 from the {N — l)-th event in X2 until the time at 
which the event N happens in Xi. This term establishes 
a bridge between the continuous-time formulation of the 
phase difference [Eq. ^] and the phase difference be- 
tween events. 

From Eq. (fTO|) . one sees that the smaller (bigger) the 
time difference jT^'— T2I is the more (the less) synchronous 
the system is, which means that the phase difference 
ki ~ ^2! gets smaller (bigger). So, it is suggestive to 
consider that the difference {r\ — r^) is linearly related 
to {Tl - T^) as 



P{Wi){Tl-n) + a{i), 



(13) 



with [3 being a constant, and a{i) brings the non-linear 
terms. 

To obtain the value of the constant (3, we imagine that 
PS is about to be lost, by a small parameter change, and 
so, — T^\ approaches j{Ti). Analogously, at this 
situation, the phase difference \r\ — has the ability 
to grow one typical cycle, i.e., (ri), and therefore, the 
term a{i) in Eq. (|13|) becomes very small and can be 
neglected. Thus, we have that l3{Wi)j{Ti) = (ri), and 
we arrive at f3 — For coherent attractors, e.g. Rossler- 
type, p is approximately 2. This result is discussed in the 
Appendixes. In Appendix , we discuss how to construct 
maps using the time events rj, and in Appendix , we 
explain how to use these maps in order to obtain that 
/3=2. 

Knowing the constant /?, we put Eq. in Eq. pTjl . 
and we have that 



N-l 



JV-1 



(14) 



Using the triangular inequality and the fact that g at 
the moment is considered to be an arbitrary constant, 
with a threshold (minimal) value, we write that 



N-l 



N-l 



\P{W^) J2 in - n)\ + I ^ a(z) + i{N)\ < g. (15) 



1=0 



i=0 



< 



Equation (fT5)) can be written as | J2iLo^ ^i}) + ^{^) 
g - \P{Wi)Y.fSQ^{Tl - T^)\. At a specific event, may 

the variable \ Yl!i=o^ '^i^) + reaches the permit- 

ted maximum value, this implies that the variable 
\[}{Wi) YhSq{TI - T^)\ gets close to zero. At this situ- 
ation, \Y^^=o ^{''')\ — Q- Using the same arguments we 
arrive at |/3(VFi) Y^iJ^ [Tl - T^)\ < g, which implies that 
(^1) < Q- Since | X^ilo^ '^(*) + max^lN)] < g we also 
have straithforward that | X]i=o^ I — Q- 

These results shows that the upper bound for the phase 
difference is given by the constant (ri) = (W) x (T). This 



means that the arbitrary constant g in Eq. ([T]) is always 
greater than or equal to (^i), in other words, (ri) is our 
threshold. The physical meaning is obvious. If (ri) is the 
bound for phase difference, given a number k > 1, the 
value K{ri) is also a bound, but it is not a minimal one. 
Thus, we fix the constant g as 



(16) 



From Eqs. 
ties 



15]) . and ((T]), we have the following inequali- 



Af-l 



i=0 

f3{W^)\Y,{Tl-n)\ < in) 



i=Q 
N-l 



i=0 



iit)~Mt)\ < in), 



(17) 

(18) 

(19) 
(20) 



If one wants to use the inequality in Eq. (|20p [or Eq. 
P9p ] to detect phase synchronization, it is required that 
the phase is an available information. For that, one needs 
to have access to a continuous measuring of at least one 
variable. The inconveniences of using this approach be- 
comes evident when either one has an experimental sys- 
tem where the only available information is a time series 
of events, like the dripping faucet experiment 17] or the 
signal is so corrupted by noise that one can really only 
measure spikes in neurons [l8| . In these two cases one 
should use the inequality in Eq. ^ [or Eq. p8| ]. The 
only inconvenience in the use of this inequality is that 
one should be careful with the type of event chosen. If 
the specified event is the spiking times, one might not see 
PS in the bursting time (and vice- versa). In detecting PS 
in large networks, it might be computationally costy to 
check for all the phase difference or event time difference 
among each pair of subsystems. In this case, one could 
check the validation of inequality in Eq. ([7]), having in 
mind that such a condition is only a necessary condition 
for PS. 



PHASE SYNCHRONIZATION IN THE TWO 
COUPLED ROSSLER OSCILLATORS 



To illustrate our approach we consider two non- 
identical coupled Rossler oscillators, given by 



2)1,2 = ai,2a;i,2 + ayi,2 

Zl,2 = fe+ ^1,2(2:1,2 - rf), 



(21) 



with ai = 1, and a2 = ai + Sa. The constants a=0.15, 
6=0.2, and d=10 are chosen such that we have a chaotic 
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attractor in a phase coherent regime. The subspace 
where the phase is computed is given by Xi = {xi,yi) 
and X2 = {x2,y2)- The time series, that define 

the events in Xj, are defined as follows. M^: tI is con- 
structed measuring the time the trajectory crosses the 
plane ?/2=0 in X2; '''2 constructed measuring 

the time the trajectory crosses the plane yi=0 in Xi. 
In Fig. [1] we show the coupled Rossler oscillator for 
the parameters 5a — 0.001 and e=0.01. We show that 
Eq. dHl) is always satisfied (for 10^ pairs of events), i.e., 
\Eto(Tl~m < (Ti)/2, with (Ti)/2=3.0353. 



0.5 <Tj> = 3.0353 




20000 40000 60000 80000 100000 

N-th event 




FIG. 2: In (A), we show the phase difference at the time 
the A''-th event happens in both subsystems. In (B) we show 
S,{N), and in (C), we show the phase difference, at the time 
that the A'^-th event happens in Xi. So, the number of events 
in X2, N2, can assume either one of the following values [N — 
\,N,N + 1). 5a = 0.001 and e=0.01. The phase is calculated 
from Eq. Q 



FIG. 1: The fluctuation | E^o(ri' - T^)], in Eq. Q. Note 
that I E^o(T'i -^2)1 < 0.5(ri). 5a = 0.001 and e=0.01. The 
phase is calculated from Eq. ([3]). 

In Fig. UK A) , we show the phase difference at the time 
the A''-th event happens in both systems, i.e. the term 
EiIo('"i - ^2) in Eq. Note that the time that the 

A^-th event happens in Xi is different that the time the 
A'-th event happens in X2. In (B), we show ^(A^) in Eq. 
P^ . and in (C), we show the phase difference, at the 
time that the A'^-th event happens in Xi. Note that the 
phase difference in (C) is just the phase difference for the 
same number of events [in (A)] plus the term ^(A^) [in 
(B)]. 

Then, we show in Fig. [3] that the linear hypothesis 
between r\-r\ and {Wi){Tl-Ti) done in Eq. ^ stands 
and /3=2.0512±0.0003. If PS is not present, such Hnear 
scale is not anymore found for the system considered. 

In Fig. ©(A), we show the quantity a in Eq. (flT)) 
for a situation that PS exists. As we decrease the cou- 
pling, Eq. (PO)) is not anymore satisfied as shown in (B), 
as well as, Eq. ([TT]). In (C) we make a zoom in of the 
vertical axis. Note the different nature of the fluctua- 
tions of the phase difference in (B) and the term a in Eq. 
p?)) . That is because the term a represents the phase 
difference without the linear growing trend, responsible 
to make the phase difference in (B) to present a positive 
slope. 

In order to compare the phase as defined in Eq. ([3]) 




<Wj> (t; - T^) 



FIG. 3: The variable r\ - 4 versus {Wi){Tl - T-l). We find 
that ri - r| ~ (3{Wi}{Tl - TJ), with /3 = 2.0512(3). 

(for Sa^O.QQl and e=0.01), and the phase as defined in 
y, e.g. 9 ~ tg{y/x), we compare the function (rj), as 
calculated for both definitions. For the phase, as defined 
in Eq. ([3]), we arrive at (rj)=6.2984 and so, (r^) > 27r. 
Other quantities are (Wj)=0.1651, and (T^) = 6.07097. 
On the other hand, the phase as defined in Q is a func- 
tion that grows in average 27r each time the trajectory 
crosses some Poincare section, which gives {vj) = 2tt. 
So, the phase definitions arrive at two different quanti- 
ties, but Eq. (PO)) is valid in order to detect PS and Eq. 
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FIG. 6: Discrete phase difference — ^2)! for no 

coupling (A) where there is not PS and for a coupUng e = 
0.00535 (B) responsible to induce PS. 



FIG. 4: In (A), we show the quantity a in Eq. (fT7|l for a 
situation when PS exists. As we decrease the coupling, Eq. 
(|20|) is not anymore satisfied as shown in (B), as well as Eq. 
(|17|) . as shown in (C). In (C) we have made a zoom in of the 
vertical axis. In (A), 5a — 0.001 and £=0.01 and in (B-C), 
Sa = 0.001 and e=0.000001. 



a parameter mismatch of Sa = 0.0002, and for a null cou- 
pling, e—0, both Rossler oscillators (presenting the fannel 
attractor) are not phase synchronized as one can check 
in Fig. [SUA), which shows the absolute discrete phase 
difference in Eq. As we introduce the coupling e = 

0.00535, the oscillators presents weak coherent motion. 



© is valid to measure the average angular frequency of 
the attractor, using both phase definitions. 



PHASE SYNCHRONIZATION IN TWO 
COUPLED NEURONS 





20 -20 -10 



FIG. 5: Chaotic attractors projected on the variables xi and 
2/1. (A) The coupling is null and therefore, there is no PS. 



Here one sees the non-coherent Fannel attractor. (B) 
coupling induces PS, creating a coherent dynamics. 



The 



To illustrate the generality of the phase definition in 
Eq. ([3]) in order to detect the phenomena of PS in non- 
coherent attractors, we consider Eqs. ((2T|) with the fol- 
lowing set of parameters, a=0.3, 6=0.4, d—7.5, such that 
we have the fannel attractor shown in Fig. ETA). This at- 
tractor has a non-coherent phase character 20, 21 1. For 




^ (A) 



5 000 5 200 5400 5600 5800 




(B) 



FIG. 7: A sample of the variables xi{n) and X2{n), from the 
subspaces that correspond to both neurons, for a situation 
where there is PS, for e=0.03 (A), and for a situation where 
there is not PS, for e=0.001 (B). In (A), Eq. (O is satisfied, 
and in (B) is not. In (A), we show three bursts, which are 
basically a sequence of spiking. 

Now, we give an example for the detection of PS with- 
out the knowledge of the state equations, but either only 
using a time series of bursting events. We consider two 
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Xj(n) 



FIG. 8: In (A-B), we show the absolute difference between the 
time of the A'^-th burst, in both subspaces (that represent the 
two neurons). In (A), Eq. © is satisfied with (Ti) =259.028 
and, in (B), Eq. © is not satisfied (there is not PS). For (A) 
and (C), Qi=4.99 and e=0.03. For (B) and (D), ai=4.99 and 
6=0.001. 

non-identical coupled neurons described by 

x,{n + l) = f[x,{n),y,{n)+(5,{n)] (22) 
yj{n + l) = y,{n)-e{xAn) + l)+ea,+ei3,{n), (23) 

which produces a chaotic attractor, for 6 — 0.001, a2=5, 
CTi = 0.240, and cr2=0.241. The subspaces are defined as 
Xj = {xj,yj). Pi^2in) = g[x2,iin) - xi^2in)]- The function 
/ is given by 

/ = «j/[l"2;j(n)]+yj(n), x<0 
f^aj+ yj{n), 0<x<aj{n)+y (24) 
/ = -1, x>a + yj{n) 

The control parameters are ai and g, with \ai —012 \ being 
the parameter mismatch and g the coupling amplitude 
(cf. The time at which events occur is defined 

by measuring the time instants in which the variable Xj , 
of subsystem Xj^ is equal to Xj ~ —1.2 (the event is 
the occurrence of a burst) [l^l , and Nj is the number of 
bursts of the subsystem X^. In this example, PS exists 
if Eq. ^ is satisfied, which also means that Eq. ([7]) is 
satisfied. 

In Fig. [71 we show the variables x\(n) and X2(n) for 
a situation where there is PS (A), and for a situation 
where there is not PS (B). Note that in (A), although the 
neurons are phase synchronized, the difference between 
the number of bursts in the variable x\(n) minus the 
number of bursts in the variable X2 (n) might be different 



than zero (for a short moment), as the hypothesis done 
in Eq. ([7]). In (A), we also represent by the dashed 
line the threshold, Xj = —1-2, from which the events are 
specified. 

In Fig. [SJA-B), we show the absolute difference be- 
tween the time of the iV-th burst, in both neurons. In 
(A), Eq. dni) is satisfied (there is PS) with (ri)=259.028, 
and therefore 0.5(Ti)=124.5014, much bigger than the 
maximum fluctuation in (A). In (B), there is no PS. In 
(C) and (D), we show a projection of the attractor on the 
variables (xi,2/i). In (B) and (D), (Ti) = 396.964, and 
(Ta) = 398.407. 

Note that although the attractor of these neurons have 
not the dynamics of a limit cycle, presenting a very com- 
plicated geometry in the phase space (as one can see in 
Fig. [8]), it is still possible to well define events as well as 
the average period of the spiking times by the use of the 
threshold shown in Fig. [71 a characteristic that defines 
this attractor to be of the weak coherent type. 

CONCLUSIONS 

We estimate the inferior bound value of the absolute 
phase difference between two coupled chaotic systems, in 
order to verify the existence of phase synchronization be- 
tween them. Our approach shows that this bound value 
(r) is given by the average evolution of the phase, calcu- 
lated in a subspace of the attractor, for a series of pairs 
of events in this same subspace. These events can be the 
number of local maxima or minima in the trajectory, the 
crossing of the trajectory to some Poincare section, or 
the occurrence of a burst/spike. 

This result was achieve because we can inspect for 
phase synchronization looking for the phase difference 
at the times for which the same number of events hap- 
pens in both subsystems. The advantage of looking at 
the phase difference at these times, instead of looking at 
the continuous phase difference, is that this approach al- 
lows us to detect phase synchronization by looking for 
a bounded time difference between events. This is help- 
ful for chaotic systems from which there is no available 
information about the state equations. 

If only the number of events is available, one can also 
find evidences of PS by checking that the absolute dif- 
ference between the number of events has to be smaller 
or equal than 1. This allows to infer the existence of PS 
in maps. These maps can be derived either from a fiow 
(by measuring events as the trajectory crosses a Poincare 
section, or by detecting local maxima of the trajectory) 
or they can be dynamical systems with a discrete formu- 
lation. 

In this work, the phase is introduced to be a quan- 
tity that measures the amount of rotation of the tangent 
vector of the flow. 

All our results are extended to coupled chaotic systems 
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that present coherent properties as defined in [20|, i.e., it 
is possible to define an average time between two events 
(T), such that for each returning time r, it is true that 
|r - (r)| < K, with K < (T). 
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CONTRUCTING PS-SETS FROM THE EVENT 
TIME SERIES 



The event time series r* can be used to construct maps 
of the attractor, whose geometrical properties states 
whether there is PS. These maps are constructed follow- 
ing simple rules: 

• At the time T2, a point of the trajectory in Xi is 
collected. 

• At the time r^J, a point of the trajectory in X2 is 
collected. 

So, as a result of measuring the trajectories in Xi 
(resp.A'2) at the times T2 (resp r^) we have a discret set 
of points T>i ( resp. 1)2) ■ 

In PS, these sets Vj will be localized, not spreading 
out to the whole attractor. In this case Vj is called PS- 
set. The theory for characthering and constructing these 
sets is presented in 13[. In a short, what happens is 
when phase synchronization occurs, the 



the following: 

times for a trajectory to pass through a Poincare section 
(or reach some defined event) becomes relatively more 
regular. Since we measure the trajectory on Xi by the 
timing of events in X2, these maps are localized around 
the Poincare section. For a non synchronous phase dy- 
namics, the sets Vj spreads over Xj. Thus, by detecting 
a PS-set, one does not have to check if the inequality for 
the time event difference holds. 

In the following, we give examples of the PS-sets 
in the coupled Rossler oscillators and in the Rulkov 
map, that mimics the neuronal dynamics presenting spik- 
ing/burting behavior. 



In Fig. [9l we show the coupled Rossler oscillators for 
a situation where PS exists. In this figure, we show bidi- 
mensional projections on the variables of subsystem X2 
(A) and Xi (B). In gray, we show the attractor projec- 
tion, and in black, projections of the PS-set T>i (A) and 
2?2 (B). Note that the PS-sets, do not visit everywhere 
Xj, rather are localized structures. 




"1 — ' — I — ' — r 



(B) 




FIG. 9: Bidimensional projection of the attractor (gray) and 
of the projections Dj of the PS-set (black) on the subspaces 
Xj. The PS-set projection ©2, in (A), is constructed using 
time series M2, and Vi, in (B), is constructed using time 
series Ml- 5a = 0.001 and e=0.01. 



PS-sets for the coupled Rulkov Map 

In the neuronal dynamics is not possible to define a 
Poincare section, due to the non-coherence of the attrac- 
tor. However, it is possible to define an event where the 
dynamics is weak coherent. This event is the ending or 
the beggining of the bursts, and in here we choose the 
beggining of the burst. Hence, we construct our time se- 
ries by measuring the crossing of the trajectory with the 
threshold, Xj = —1.2. 

In Fig. [iniwe show a projection of the attractor on the 
variables {xi,yi), in black points, and the subsets Vj, in 
gray points. In (A), where we have phase synchronization 
the set Vi does not fuUfil the whole attractor, but is 
rather localized, whereas in (B), where PS is not present 
the set I?2 spreads over the attractor X2. 



PS-sets for the coupled Rossler oscillators 

The time series, A^J, that define the events in Xj, are 
defined as follows. Ml: t{ is constructed measuring the 
time the trajectory crosses the plane ?/2=0 in X2; M^: 
is constructed measuring the time the trajectory crosses 
the plane yi=0 in Xi. 



P DIGRESSION 

In this section we explain why in coherent attractors, 
e.g. Rossler-type, the constant /? is approximately 2. 

That is so, because we compare the phase difference at 
the time events occurrence. Let us just remember that 
we are measuring the phase in one subsystem at the times 
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FIG. 10: In (A-B), we show the absolute difference between 
the time ol the A'^-th burst, in both subspaces (that represent 
the two neurons). In (A), Eq. © is satisfied (there is PS), 
with (Ti) =259.028 and, in (B), Eq. Q is not satisfied (there 
is not PS). For the same parameters as in (A), in (C) we show 
the PS-set projection I>i (gray), while for the same parame- 
ters as in (B), in (D) there is not a PS-set. For (A) and (C), 
ai=4.99 and e=0.03. For (B) and (D), qi=4.99 and e=0.001. 



(A) 




(B) 




FIG. 11: Pictorial visualization of a situation where PS exists. 
We represent in (A) the trajectory in the subspace Xi, and 
(B) the trajectory in the subspace X2. The filled regions in 
(A) [(B)] represent trajectory positions at the time the Ai'-th 
event happens in X2 [Xi] for a situation when PS exists. An 
event is considered to happen when the trajectory crosses the 
dotted line. 



that events in the other subsystem happen. Hence, at the 
time events happens in Xi [resp. A^j], we collect points 
in X2 [resp. Xi], obtaining the gray filled region in Fig. 
fTlTB) [resp. (A)], which is a PS-set. 

In particular, when the A'^-th event happens in X2, the 
trajectory on Xi is indicated by the cross in (A) . At this 
time, , we record the phase in Xi, namely ?!'i(t^). 
As the time goes on, the trajectory (in a unclockwise 
direction of rotation) on Xi reaches the event line in X2 
at the time . At this time, the trajectory in X2 is at 
the cross in (B), and the phase is (I)2{t^)- Since these are 



the particular case represented in this figure. That is so 
because the time difference is approaximately given by 
the time that the trajectory in Xi spends from the cross 
in (A) till the event line, which is approximately 1/4 of 
the average period (T) . 

The phase difference, at which the same number N of 
events happen, is [0i(rj'^) — (piir^)] ~ (^)/2, since this 
phase difference is basically given by the displacement 
of the phase in Xi from the cross in (A) till the event 
line, plus, the displacement of the phase in X2 from the 
event line till the cross in (B). But that is approximately 
given by 1/2 of the average increasing of the phase (r), 
which was shown to be equal to (W) x (T). Therefore, 
\4'i{tI') ~ (j>i{T^)\ « 2 X \taU2 —t(^\, which consequently 
give us that « 2. 



typical events, we can say that [tj — t 
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